Dynamical renormalization group approach to the Altarelh-Parisi-Lipatov equations 
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The Altarelli-Parisi-Lipatov equations for the parton distribution functions are rederived using 
the dynamical renormalization group approach to quantum kinetics. This method systematically 
treats the InQ'^ corrections that arises in perturbation theory as a renormalization of the parton 
distribution function and unambiguously indicates that the strong coupling must be allowed to run 
with the scale in the evolution kernel. To leading logarithmic accuracy the evolution equation is 
Markovian and the logarithmic divergences in the perturbative expansion are identified with the 
secular divergences (terms that grow in time) that emerge in a perturbative treatment of the kinetic 
equations in nonequilibrium systems. The resummation of the leading logarithms by the Altarelli- 
Parisi-Lipatov equation is thus similar to the resummation of the leading secular terms by the 
Boltzmann kinetic equation. 

PACS numbers: ll.lO.Hi, 12.38.Aw, ll.lO.-z 

An early important test of QCD and one that validates the parton picture, is the evolution of the quark 
distribution, or structure function, in deep inelastic lepton-hadron scattering which reveals scaling violations as 
predicted by QCD. The leading logarithmic dependence of the quark structure function on the scale in deeply 
inelastic scattering was originally computed in QCD via the operator product expansion |l| . In perturbation theory, 
parton emission leads to QCD corrections to the distribution functions that grow as In Q , which eventually become 
large and invalidate the perturbative expansion. These corrections signal the violations of scaling in deep inelastic 
scattering. 

Lipatov and Altarelli and Parisi Q provided an alternative derivation of the evolution equations for the parton 
distribution functions by associating the In corrections with a change in the parton distribution function at the 
scale . These authors provided a set of integro-differential equations that resum the perturbative expansion and 
describe the evolution of the parton densities as a function of In Q . These parton evolution equations are successful 
in predicting the correct leading logarithmic dependence of the distribution functions in and the scaling violations 
predicted by these equations have been confirmed spectacularly in deep inelastic experiments . 

Altarelli and Parisi used a probabilistic description in their derivation, however infrared singularities associated 
with the emission of low momentum partons introduce singularities in some probability densities. Anticipating that 
the infrared singularities should cancel between real and virtual soft gluons, these authors provided a prescription to 
handle these singularities. 

Collins and Qiu |^ provided a derivation of the Altarelli-Parisi-Lipatov (APL) equations with a special treatment 
of virtual diagrams that clearly reveals the explicit cancellation of the infrared singularities between the contributions 
of real and virtual soft gluons through Ward identities. The result of Collins and Qiu is contained in a set of integro- 
differential equations, equivalent to the APL ones but without the need of extra prescriptions for handling soft gluon 
emission since the cancellation of infrared singularities is manifest in their treatment of real and virtual diagrams. 
For valence quarks these equations have the form 
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where qi{x, Q^) is the distribution function of the (valence) i-quark at a fixed momentum fraction x that evolves with 
the momentum scale at which the distribution function is probed and the quark-quark splitting function is given 

by 
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with Nc the number of colors. The splitting function ^qq{z) describe the probability density of finding a parton 
quark inside another parton quark with momentum fraction z of the parent parton momentum and was originally 
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introduced by Weiszacker and Williams for QED 0| . The relative minus sign between the two terms allows an explicit 
cancellation of the infrared divergences arising from (real and virtual) soft gluon emission |^ . 

A remarkable aspect of the Altarelli-Parisi-Lipatov equations in the form obtained by Collins and Qiu given by 
Eq. (0) above, is that they feature exactly the same splitting function in both terms on the right-hand side. This 
feature allows an unambiguous probabilistic interpretation as a set of Boltzmann kinetic equations^ since the different 
contributions can be interpreted as "gain" and "loss" processes that change the parton distribution in "cells" labelled 
by the variable x, the fractional momentum of the hadron carried by the parton 

Within the context of the probabilistic interpretation suggested by Collins and Qiu and for reasons that will become 
clear below, it is important to remark at this stage that the Altarelli-Parisi-Lipatov equation involves a convolution 
in the variable x but not in the variable Q^, i.e., the parton distribution function on the right-hand side of the 
equation is evaluated at the same as that on the left-hand side. This feature of the evolution equation results in 
a Markovian approximation of parton evolution inside a hadron. Physically, this implies that the evolution of the 
parton distribution at some fixed x involves only parton emission or absorption processes taking place at momentum 
scale where the parton distribution is probed, rather than those taking place at any other momentum scales. This 
property is actually inherent from the leading logarithmic approximation employed in the formulation, where each 
parton emission or absorption process is actually independent of the others and depends only on the momentum scale 
at which the process takes place. 

This interpretation of the Altarelli-Parisi-Lipatov equations as a set of kinetic equations for the parton distribution 
functions motivates us to provide yet another alternative derivation of these equations by implementing a novel 
approach to quantum kinetics introduced recently in the form of the dynamical renormalization group 

While we anticipate that the approach presented here will ultimately lead to the form of the Altarelli-Parisi-Lipatov 
equations derived by Collins and Qiu, we believe that this alternative formulation yields novel insights that could prove 
valuable to implement a systematic treatment of parton evolution. As will be shown in detail below, the rederivation 
of the evolution equations via the dynamical renormalization group provides further and deeper understanding of 
many important aspects: 

(i) It reveals some subtle aspects of the derivation which highlight the role of the running of the strong coupling 
and the "separation of scales" which is a necessary ingredient in a kinetic description. The role of the running 
coupling is ultimately at the heart of the scaling violations and the necessity for including the running of the 
coupling in the evolution kernel becomes manifest in this approach. 

(ii) The dynamical renormalization group directly reveals how the InQ^ that arise in the perturbative expansion 
can be systematically absorbed into the parton distribution functions and establishes a direct correspondence 
with secular divergences in time Q that are manifest in the derivation of the Boltzmann equation and also 
responsible for the failure of the perturbative expansion in quantum kinetics. 

(iii) It directly establishes a correspondence between the leading logarithmic approximation to parton evolution 
equation which is local in and the Markovian approximation to kinetic equation which is local in time. 

As noticed by Collins and Qiu and mentioned above, the connection between the Altarelli-Parisi-Lipatov equation 
and the kinetic equation in nonequilibrium statistical physics is established by identifying In with the "time" 
variable. This identification is the main concept upon which the formulation based on the dynamical renormalization 
group hinges. In this article, we revisit the derivation of Eq. but implementing the concepts of the dynamical 
renormalization group |9| . 

The i-quark distribution function at some scale with fixed momentum fraction x is defined as Q 

q.{x,Q')^i2n)'2E,l ^l-TUP,?) {l-)ba S [x ~ ^) , (3) 

with 

TUP,p) = J^J d^ye-^P-y{P\My)MO)\P), (4) 

where = (n+, n_, n±) = (1, 0, 0_l) and P is the momentum of the hadron. The derivation of the parton evolution 
equations to a given order in the strong coupling, begins by expanding T^^ in perturbation theory in as, and evaluating 
the corresponding Feynman diagrams at large and P with x — Q^/2P ■ q fixed. 

The main idea followed in Ref. |^ is to calculate the change in the parton distribution function from ghion radiation 
processes. In the infinite momentum frame as the scale changes the number of partons with momentum fraction x 
also changes. Consider the number of partons with momentum fraction x at the scale Q^. On the one hand, partons 
with larger momentum fraction xi > x a.t the scale that emit a gluon with momentum fraction xi — x now enter 



3 



into the "phase space cell" with fractional momentum x. On the other hand, partons that are already in the "phase 
space cell" with momentum fraction x at the scale leave this "cell" by diminishing their value of x upon radiating a 
gluon. The former corresponds to the "gain" contribution, whereas the latter corresponds to the "loss" contribution. 

Collins and Qiu showed that the real gluon emission processes gives a positive contribution to the parton distribution 
function and can be calculated from the real gluon diagram by introducing a cutoff Q. Most importantly, by 
converting the virtual gluon diagrams into a real diagram which has the same form as that corresponding to real gluon 
emission except for having the opposite sign, they obtained an explicit cancelation of infrared divergences between 
diagrams and the "gain" and "loss" contributions to the parton distribution function 



Agi(a;, A^) ^ Aqi{x, A^)g^in - A9,(a;, A2)ioss, (5) 
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where the quark-quark splitting function ^qq{z) is given by Eq. (|^). Taking the derivative with respect to InA^ and 
identifying the cutoff scale A^ as the scale where the parton distribution function is probed , Collins and Qiu 
obtained the parton evolution equation for valence quarks as given by Eq. (|l|). 

A close examination the above equations reveal the following important features in their derivation: 

(i) The quark distribution function that enters the right-hand side of the gain and loss contributions given by 
Eq. carries the k\ dependence. Since k\ is integrated out to give the total change A(7i(a;, A^), this in turn 
is tantamount to an implicit resummation of the leading logarithms in Collins and Qiu's derivation (see also 
below) . 

(ii) The strong coupling is taken to be running with the scale of the exchanged gluon in the process. While the 
lowest order calculation of the relevant Feynman diagram given by Fig. 5 in Ref. ^ would be in terms of the 
bare gluon propagator and the bare strong coupling constant, in allowing the strong coupling to run with the 
scale, the exchanged gluon has been effectively "dressed" by parton loops. 

We address these issues in detail by implementing the novel dynamical renormalization group approach to quantum 
kinetics This approach reveals in an alternative manner, the type of resummations implied by the final form of the 
parton evolution equation. This alternative formulation is based on Collins and Qiu's probabilistic interpretation of 
the parton evolution equation as a Boltzmann kinetic equation with In interpreted as a "time" variable and a; as a 
"momentum" or "cell" variable in phase space. In particular, as will be shown below, the dynamical renormalization 
group approach to the Altarelli-Parisi-Lipatov equations not only provides and justifies an alternative and explicit 
resummation of the leading logarithms in the form a resummation of the leading secular terms, but also indicates 
unambiguously the running of the strong coupling in Collins and Qiu's derivation is crucial to violations of scaling. 

The derivation of the kinetic equation for the single-particle distribution functions in phase space via the dynamical 
renormalization group |^ begins by considering that the distribution function is determined at some initial time 
to at which the initial state of the nonequilibrium system is specified. That is, the initial condition for the time 
evolution is determined by the initial distribution function nk(to) for particles of momentum k. Time evolution 
through the interaction, changes the number of particles in a given phase space cell, say, labeled by momentum k. 
The evolved distribution function at a later time t = <o + At is calculated in perturbation theory and can be written 
as nk(i) = f^k(io) + Ank[nk(io), At]. The change in the distribution function depends on the distribution function 
at the initial time nk(to) because the perturbative expansion is carried out with respect to the initial state the 
dependence on At is determined by the transition probability which is obtained by using time-dependent perturbation 
theory. For At much larger than any microscopic time scales, using Fermi's golden rule one finds that the transition 
probability is proportional to At so that Ank[nk(to), At] = C[nk(to)] At, where the collision term C[nk(to)] depends 
on the distribution functions at tg. Hence, one obtains 

^=C[.k(to)], (7) 



which is almost the form of the Boltzmann equation, were it not for the fact that the term on the right-hand side 
depends on the distribution function at the initial time. This in turn implies that the change in the distribution 
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function at long times is proportional to At. Thus perturbation theory features secular terms, i.e., terms that grow 
in time and would lead to the breakdown of the perturbative expansion. 

The Boltzmann equation is obtained after a resummation directly in real time via the dynamical renormalization 
group The resummation implied by the dynamical renormalization group can be summarized as follows: Consider 
evolving in time from to to to + At with At much larger than a typical microscopic time so that Fermi's golden 
rule can be applied, but not too long so that perturbation theory is still valid during this time interval. This is 
exactly the assumption of the separation of time scales in a kinetic description. At the end of this interval the new 
distribution function is calculated as nk(to + At) = nk(to) + C[nk(to)] At. Now use this new distribution function 
as the initial value at to + At to evolve further another time step, i.e., the distribution function is "reset" to the 
new value after each evolution in a time step. This resummation procedure can be implemented consistently via the 
dynamical renormalization group |^ and leads to the Boltzmann kinetic equation 

jn^{t) = C[n^it)], (8) 

whose solution is found to be given by 

nk(t) =nk(to) + Ank(t), Afii^it)^ f dt' C^')]. (9) 



A comparison of Eq. (||) and the gain and loss contributions to the change of the parton distribution given by Eq. (H) 
reveals clearly that a resummation is implicitly invoked in Collins and Qiu's derivation. 

We now implement a resummation via the dynamical renormalization group to obtain the evolution equation for 
the parton distribution function in the same manner. Following Ref. ^ we now focus on the evolution equations for 
the distribution function of valence quarks resulting from gluon emission. Those for the sea quarks can be obtained 
by a straightforward generalization once the relevant Feynman diagrams have been computed. 

Consider that the parton distribution qi{x, Qq) at a scale Q§ is given and that the evolution from the scale Qq to a 
larger scale results in a change of the parton distribution as a consequence of parton radiation. These changes are 
calculated via a perturbative expansion in the strong coupling constant by following exactly the steps summarized 
above and described in detail in Ref. |^. The freedom to choose the factorization and renormalization scale permits 
to identify the upper momentum cutoff A'^ as in the loop integrals leading to Aqi{x,Q^) B. Taking the lower 
limit to be the scale Qo ^^'^ distribution function that enters in the diagrams at the scale consistent with the 
leading logarithmic approximation, we obtain 

qi{x,Q^) = qi{x,Ql) + Aqi{x,Q^;Ql), 

Aq,ix, Q2; Ql) = Aq,{x, Q2; Qg)g^i„ - Aq,{x, Q^; g^)ioss, (10) 



where 



Agj(x,(5^;Qo)gain = / TrTr"''('^i) / Qli^l, Ql) Iqq , 

Agj(a;,(5'^;Qo)ioss = / 7r^ois{k\)qi{x,Ql) i dz-fqq{z), (11) 
j2 znn^ Jq 
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"^(^-) = ln(fci/A^,,)' "°-3332]v7' 



and Nf is the number of flavors. 

We emphasize an important aspect of the above equations that is different from those obtained by Collins and Qiu: 
the parton distribution function enters at the initial scale Qq. Thus in the "update" equation (|l^) with the gain 
and loss contributions given by Eq. ( pi] ) the parton distribution function is only convoluted in the variable x but not 
in fc^. This aspect in turn allows to obtain the characteristic scales of parton evolution already at the perturbative 
level. The dependence of the updated parton evolution on the initial parton distribution 9i(a;, Qg) is reminiscent of 
the derivation of the Boltzmann kinetic equation from the time evolution of a given initial state of a nonequilibrium 
system outlined above. 

Although, following Collins and Qiu, we have included the running of the strong coupling in the gain and loss terms, 
we now study in detail the justification for this choice. Consider first the above expression in strict perturbation theory 
by taking the strong coupling to be a constant q;s(Qo) fixed at the initial scale Qq. This corresponds to considering 
a bare exchanged gluon in the corresponding Feynman diagram. Performing the k'^ integral, one obtains 

Aq,{x, Ql) cx asiQl) HQVQD- (12) 
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In the interpretation of In as the "time" variable t, this in turn translates into the statement 

Agj(a;, i; t^) oc ^^(^o) - ^o) , (13) 

which is akin to the expression for the real-time evolution of a single-particle distribution function calculated in strict 
perturbation theory by using Fermi's golden rule as described above. 

In the dynamical renormalization group approach |ll[ these type of terms (divergences in the long-time limit) 
are identified as secular divergences in the perturbative expansion which signal the breakdown of perturbation theory 
and can be resummed systematically via the dynamical renormalization group. 

Equation (|l^) reveals that the parton distribution function evolves on "time" scales (not to be confused with the 
real time) comparable to the scale of the running of the strong coupling constant. Therefore during the "time" scale 
of evolution of the parton distribution the coupling actually changes substantially and its "time" evolution cannot be 
neglected. This is to say that keeping the strong coupling as a constant (determined at the initial scale) does not allow 
a clear separation of "time" scales between the "kinetic" scales characteristic to parton evolution and the "microscopic" 
scale determined by the strong coupling since they are of the same order in the logarithmic approximation. 

Including the running of the strong coupling in the loop momentum integral is tantamount to including the vacuum 
polarization in the exchanged gluon and leads to 



A(ji(a;, Q^; Qo) °^ In 

or, in terms of the "time" variable t, 



(14) 



Agi(a:, i; to) oc ln(t/to). (15) 

Therefore, including the running of the coupling in the evolution kernel leads to a slower evolution of the parton 
distribution function. This is clearly a consequence of the fact that during the "time" scale of evolution of the 
distribution function the strong coupling changes substantially, becoming smaller and therefore in turn slowing down 
the change in the parton distribution. 

This is the important issue of "time" scales that emerges in the kinetic description. If the coupling varies on "time" 
scales longer than that of parton evolution, then it can be kept constant during the evolution since there is a wide 
separation of "time" scales. However, if the coupling varies substantially during the "time" scale of parton evolution, 
then the running of the coupling must be included in the "update" equation of the parton distribution function. 
Hence, this formulation based on the "update" equation (|l^) with the gain and loss contributions given by Eq. ( |Tl| ) 
in terms of the parton distribution function at the initial scales allows to clarify unambiguously whether or not to 
include the running of the strong coupling by comparing the "time" scale of parton evolution and that of the running 
coupling. The criterion is as follows: If the scales are comparable or the parton distribution evolves slower ^ then the 
running of the coupling must be included; if the parton evolution evolves much faster, then the coupling can be kept 
fixed. 

In the case of evolution of the parton distribution functions in QCD, the running of the strong coupling must be 
included in the evolution kernel of the "update" equation. A similar situation in which the evolution kernel depends 
explicitly on time has been found in studying quantum kinetic of charged particles in a hot QED (QCD) plasma 
as a result of infrared divergences associated with the emission and absorption of quasistatic, soft magnetic photons 
(gluons) [|| The similarity between these cases is tantalizing. 

Furthermore, the above "update" equation (|o|) with the gain and loss contributions given by Eq. (|l^) implies that 
for fixed x the change in the parton distribution function due to a change in momentum scale only depends on the 
parton distribution at the initial scale Qq, i-e., the momentum scale at which the gluon is emitted, rather than on 
the parton distribution at intermediate scales between Qq and . We want to stress once again that this result 
is consistent with the leading logarithmic approximation where gluon emission at each vertex is independent of the 
others and depends only on the momentum scale at which the gluon is emitted. It is this feature of the "update" 
equation that establishes a correspondence between the leading logarithmic approximation in parton evolution and 
the Markovian approximation in a kinetic description ||ic| . 

The logarithmic divergences for large either in the form of Eq. (^2|) or (|l^) clearly signals the breakdown of the 
perturbative expansion and the necessity for resummation. These divergences are akin to the secular divergences that 
appear in time-dependent perturbation theory, in the derivation of kinetic equations ||^ , as well as in the asymptotic 
analysis of differential equations |pT| . A systematic method based on the idea of renormalization group has been 
developed in Ref. |ll| to obtain asymptotically convergent solutions of differential equations whose perturbative solution 
features secular terms. This renormalization group method, which we now follow, has been recently generalized to 
the realm of nonequilibrium quantum field theory and quantum kinetics [p[ . 
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The main idea of the renormahzation group method, as described above in obtaining the Boltzmann equation, 
is to evolve the parton distribution function from an initial scale Qq in perturbation theory up to a larger scale 
Qq + AQ^. The change in scale AQ^ is small enough so that the change in the parton distribution is perturbatively 
small and hence the perturbative expansion is reliable. The change in the parton distribution is absorbed into the 
updated parton distribution at the scale Qq + AQ^, which is then used for the next evolution in scale, and so on. This 
procedure is implemented systematically by introducing a "renormahzation constant" of the initial parton distribution 
Z(x,t'^,Qq) at an intermediate but otherwise arbitrary scale t^, i.e., 

qt{x, Ql) = Z{x, T^, Ql) q,{x, r^), 
r^, 0^) - 1 + 7S^\x, r^ Ql) + z^^\x, r^ Qg) + • ■ • 

where (7^(0;, r^) is the "renormalized" parton distribution at scale and z*^""^^ (x, r^, Qq) ~ 0{as), z^^\x,t^, 
©(a^), etc. Substituting Eq. (|l6|) into Eq. ( |Tc| ) and keeping terms only up to 0{as), one obtains 



(16) 



Jq2 27rfc^ xi 



dxi 



as{k\)qi{x,T'^] 



Ql 







dz-iqq{z) + 0{al). 



(17) 



The coefhcients z^^\x , ^ Q"^) are chosen to cancel exactly the secular (logarithmic) terms at the scale r order by 
order in the perturbative expansion. Thus, to ©(a^), one obtains 



zW(x,r2,Qg) = - 



1 



dp, 



dk\ 



as{k] 



^ dxi 

Xi 



gj(a;i,T^)7,g 



27rfc2 



as(^i)'?i(a;,r^) / dz-/gg{z) 



(18) 



Upon substituting this choice for z^^'' (x , t'^ , Qq) back into Eq. (p7|), we arrive at 



qi{x,Q^) = qt{x,T^) + 



dP, 



27rfc2 



f'^' dkl n2. dx^ , 2^ ( x\ 
^2 27rfcj^ xi \xi J 



as{k]_)qi{x,T^) dzjqq{z) + 0(0^^), 



(19) 



which can be improved by choosing the arbitrary renormahzation scale close to the scale . This renormahzation 
procedure embodies the original idea behind the Altarelli-Parisi-Lipatov equation that the change in the parton 
distribution due to a change in the scale can be absorbed into a redefinition, indeed a renormahzation, of the parton 
distribution. 

Since the scale is arbitrary, a change in the renormahzation scale is compensated by the corresponding change in 
the "renormalized" parton distribution so that qi{x, Q^) is independent of the renormahzation scale t^. This in turn 
leads to the dynamical renormalization group equation t'^ dqi{x , Q'^) / dr^ — 0, which to lowest order in is given by 



2 " / 2\ 



«s(r^) 
27r 



I —q,{xi,T'^)jqq ( — ) 
Jx \xij 



dz ^qq(z) 



(20) 



Upon identifying the renormahzation scale as the scale , one finds immediately that Eq. (^) is exactly the 
Altarelli-Parisi-Lipatov equation derived by Collins and Qiu. Consequently, we find that while the dynamical renor- 
mahzation group leads to a kinetic equation that rcsums the leading secular terms in the perturbative expansion in 
nonequilibrium quantum statistical mechanics, the Altarelli-Parisi-Lipatov equation resums the leading logarithms, 
which as argued above are identified with the secular terms in a kinetic description. 

The solutions of Eq. ( |20| ) can be found by taking moments of the parton distribution and quark-quark splitting 
functions [|| 



J^"(Q2) ^ [ dx 
Jo 



dz (z" 



(21) 



7 



After some algebra, one finds 

dQ2 ' ' 27r 



Q'4r.FTm = -^4?^7" (22) 



where 

iV2 



7, 



99 



2A^. 



1 1 



" 1 



2 n(n+l) 



(23) 



With the running strong coupling given by Eq. ( |11| ) , the solution is found to be given by 

-"07J,/2t 

m I f t- 1 1 



In(QVA^CD) 



ln(Q2/A|^^) 



(24) 



The mild violation of scaling, a logarithm raised to a power, instead of a pure power law with an anomalous dimension 
is a consequence of the running of the strong coupling inside the evolution kernel. Within the kinetic interpretation of 
parton evolution afforded by the dynamical renormalization group, this is a consequence of the fact recognized above 
at the perturbative level that the coupling runs on scales shorter than that of the distribution function and hence the 
full running coupling must be included in the kernel [see the discussion after Eq. ([Tl])]. 

In conclusion, we have given an alternative derivation of the Altarelli-Parisi-Lipatov equation from the point of view 
of the dynamical renormalization group. While the final equation of evolution for the parton distribution function is 
the same as that obtained by Collins and Qiu we believe that the derivation via the dynamical renormalization 
group provides further insight into the nature of the approximation and resummation. This approach helps clarify 
several important issues in the original derivation of Collins and Qiu: 

(i) At the level of perturbation theory, a detailed analysis of the "update" equation ( p^ ) with the gain and loss 
contributions given by Eq. (^l]) reveals that a consistent treatment of the parton evolution to leading logarithmic 
order should only involve the parton distribution at the initial scale. In turn this permits to understand 
unambiguously whether the strong coupling should be allowed to run during the scale evolution. This analysis 
is similar to the issue of separation of time scales in a kinetic description of nonequilibrium phenomena. It is 
precisely this analysis that ultimately leads to the correct violations of scaling , as a power of logarithms in 
Eq. (p^, a consequence of the running of the coupling in the evolution kernel. 

(ii) The dynamical renormalization group manifestly and systematically treats the InQ^ corrections that arises in 
the perturbative expansion as a renormalization of the parton distribution functions and establishes a direct 
correspondence between the leading logarithmic divergences arising from soft gluon emission and the leading 
secular divergences in time in the derivation of the Boltzmann equation. A common feature of such divergences 
in the parton evolution equations and in quantum kinetics is that they signal the breakdown of the perturbative 
approach to the evolution dynamics for large-momentum or long-time scales respectively. This breakdown of 
the perturbative expansion requires a resummations which is consistently recognized as a renormalization of the 
distribution functions. 

(iii) The dynamical renormalization group approach reveals that the leading logarithmic approximation in which 
the parton evolution equation only depends on the parton distribution at the momentum scale at which it is 
probed, corresponds to the Markovian approximation in a kinetic description. 

We hope that the experience with non-Markovian kinetics in nonequilibrium statistical physics can pave the way, 
via the dynamical renormalization group, to understand how to consistently treat next-to or next-next-to leading 
logarithmic corrections. Beyond the leading logarithmic approximation lack of angular ordering of gluon emission will 
likely result in convolutions both in the variable x and in the variable Q^, i.e., a non-Markovian parton evolution in 
Q2 within the kinetic interpretation. Work on these aspects is currently underway. 
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